A flux tube model for glueballs 
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We calculate the mass spectrum and the decay widths of glueballs in the flux tube model. The 
glueball is assumed to be a closed flux tube. The breathing motion and the rotational motion are 
investigated using the WKB approximation. The calculated spectra is consistent with those by 
lattice QCD. The decay widths are also computed using Schwinger mechanism and it is shown that 
they have rather large values. 



I. INTRODUCTION 



It is widely believed that the strong interaction is de- 
scribed by quantum chromodynamics (QCD) [J]. It is, 
however, much difficult to solve QCD so that several ef- 
fective models arc devised to explain physical properties 
of many hadrons. The flux tube model |2(, or hadron 
string model 0-0, is one of such models. According 
to the model, mesons are made up of a quark and an 
anti-quark connected by color flux tube and baryons are 
described by a quark and a diquark connected by the 
same flux tube. We showed in a previous paper that 
various excited states of the hadrons are explained sys- 
tematically using such a simple model. On the other 
hand, the recent lattice QCD simulations fj predict the 
existence of glueballs, which are composed of only gluons 
. Unfortunately such a hadron has not been identified yet 
experimentally but the existence seems to be probable in 
Nature. 

In this paper, we present an intuitive and analytic 
approach to the glueballs using a flux tube model. It 
is assumed that the glueball is a closed color flux tube 
(flux tube ring). Moreover the flux tube that constitutes 
the ring is assumed to have the same properties as the 
flux tube that constitutes a meson or a baryon. Con- 
sequently our model has no free parameter. We study 
the relativistic motion of the closed string and obtain 
the mass spectra of glueballs. Their decay widths are 
also calculated by taking into account qq pair produc- 
tion inside the flux tube (Schwinger mechanism). This 
simple picture has been considered by several authors. 
In particular, Koma, Suganuma and Toki 7] studied rel- 
ativistic breathing motion analytically with the use of 
the dual Ginzburg-Landau theory Q-[[3- ^ ls > however, 
very interesting to investigate the rotational motion of 
the closed string, which corresponds to the Regge trajec- 
tory of glueballs. 

This paper is organized as follows. In the next sec- 
tion, the closed flux tube is formulated and we obtain a 
Lagrangian of the string. In Sec. 3, canonical quantiza- 
tion is applied to the dynamical system and we obtain 
the Schrodinger equation. Then the mass spectra of the 
glueballs are obtained with the help of the WKB method. 



The decay widths of the glueballs are calculated by tak- 
ing into account the quark pair creation mechnism in 
Sec. 4. The numerical results are presented in Sec. 5. The 
concluding remarks are given in the final section. 



II. A STRING MODEL OF GLUEBALLS 

According to our flux tube model [5j, mesons (baryons) 
are composed of a quark and an anti-quark (diquark) con- 
nected by a color flux tube: mesons: (q — ► q) and baryons: 
(q — » qq). Each arrow denotes a color flux tube. The 
color quantum number of q is the same as that of qq 
so that both flux tubes are identical. This fact means 
the universality of the slope of the Regge trajectories of 
mesons and baryons. Of course, since the mass of q is 
different from that of qq, the lengths of the flux tubes 
are different. 

Next let us consider a closed flux tube of glueballs, 
which are supposed to be made up of two gluons. First, 
we notice that the color quantum number of a gluon is 
the same as that of a quark pair (qq) belonging to an 
octet representation of SU(3). Generally a qq pair is de- 
composed into two representations as follows: 



QQ = (qq)i © (qq)s 



(1) 



The first term in the right-hand side denotes a singlet 
representation which corresponds to a meson state. On 
the other hand, the second term which is an octet has 
the same color quantum number as a gluon. Noting that 
a quark (anti-quark) is a source (sink) of color flux, the 
gluon has two flux tubes; one of them flows out and the 
other in. Therefore if two-gluon state is a color singlet, 
two flux tubes of one gluon have to connect with other 
two flux tubes of the other gluon to form a circle (ring) 
[? ]. This closed string becomes nothing but a glueball. 
From this formation of the color flux, it is apparent that 
the flux tube of the glueball is the same as those of mesons 
and baryons. 

Now we will construct a mathematical expression for 
this glueball. Let us start with the MIT Lagrangian ap- 
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FIG. 1: The closed string with radius r rotating around :r-axis 
with angular velocity uj. 

plied to a torus (bag) [Tlj . 

L = I J(U 2 - E 2 )d 3 x - J Bd 3 x, (2) 

where E and H are color electric and magnetic fields 
respectively and the B in the second term is the bag 
constant. We consider the breathing motion of the radius 
of a flux ring and the rotational one arround a diameter 
of the ring; the closed string with varying radius r(t) 
rotates around x axis with angular velocity co = 6 (t) 
as shown in Fig.l. The volume integral in the right- 
hand side of Eq.(2) can be rewritten in the coordinates 
system which rotates together with the ring. If the color 
flux in the rotating coordinate is denoted by Eo, we get 
H 2 — E 2 = — Eq. Then the Lagrangian is approximated 
by the line integral, 

L{r,r,cj) = <j) ElASoVl - v 2 dl- <j) BASoVl ~ v 2 dl 

.(3) 

where the volume element of the ring with the infinites- 
mal length dl is approximated by ASoVl v 2 dl. The 
A5*o means cross section of the ring in the rotational co- 
ordinates system. The velocity v of the line element is 
perpendicular to the direction of the line element so that 
the factor Vl — v 2 represents the Lorenz contraction. If 
the generalized coordinates of the line element is writ- 
ten as those in the Fig.l, the squared velocity is given 



by v = 



+ r 2 u) 2 sin 2 tp. 



Consequently we obtain the 



following Lagrangian of the closed ring, 

L(r,r,ui) = — / ar\J 1 — v 2 dp, 



(4) 



where the constant a means string tension given by 
a = (E 2 /2 + B)AS . Moreover if we think of the pres- 
sure balance on the flux tube, we get Eq/2 = B 12]. The 



final expression of our Lagrangian is simple. This result 
is realized by the following: if the closed string is divided 
into many particles (line elements) Am*, the total La- 
grangian of this system becomes L = — J2i Am^ y/1 — v 2 
and this is our Lagrangian obtained above. 



III. MASS SPECTRA OF GLUEBALLS 

In this section, we quantize the Lagrangian in order 
to calculate the energy eigenvalues of the glueballs. The 
generalized coordinates of our system are the radius r of 
the flux ring and the rotational angle 9 of the ring around 
the x-axis. The conjugate momenta of these coordinates 
are defined by 



Po 



dL 
df 
dL 



Vl — V* 

or u> sin ip 



-.dip = Hi 



dp — Ilu 



(5) 
(6) 



where p r and pg represent the radial momentum and the 
angular momentum of the system. It will be also shown 
that the quantities H and I introduced in the right hand 
side correspond to the energy and moment of inertia of 
the closed string respectively. The Hamiltonian is given 
by 



H = rp r + ujpg - L 



VT 



-d(p 



(7) 



This is the energy of our system (the ring) and already 
appeared in Eq.(5). 

Now let us go over to quantum theory. The canonical 
quantization leads to replacement of the momenta by the 
following operators: 



P 



■n 1 d 
-ih-—r 

r or 
n d 



(8) 
(9) 



where the form of p r should be regarded as the three- 
dimensional radial coordinate. From this expression, the 
normalization of the wave function is 



d.r 



dipr 2 \il)(r, i 



= 1 



(10) 



Thus we obtain a Schrodinger equation with an energy 
eigenvalue E, 



H(r,p,l)ip(r,6) = Eijj(r,9) 



(11) 



Noting that [H, I] = 0, the wave function can be written 
as 



i 



27T 



ip(r)e 



ne/h 



(12) 
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Here I is an eigenvalue of the angular momentum, which 
is an even integer (I = 0, 2, 4, • • •) because of the bound- 
ary condition: ip(r,6 = 0) = ip(r,9 = 2ir). Then the 
Schrodinger equation is reduced to the equation for the 
new function tp(r) of r, 



H(r,p, l)ip(r) = E(p(r) 



(13) 



To solve this equation, let us divide two cases: I = and 
I ^ 0. 



In this case, we get lo = from Eq.(6). The Hamil- 
(2nar) 2 so that the Schrodinger 



tonian becomes 
equation is reduced to 



{p 2 + {2irar) 2 }ip(r) = E 2 ip(r) 



(14) 



It should be noted that this equation is nothing but 
Schrodinger equation of the three-dimensional harmonic 
oscillator. Introducing a new function x( r ) by = 
x(r)/r, we have the Schrodinger equation for the one- 
dimensional harmonic oscillator so that we obtain the 
familiar eigenvalues for E 2 , 



E 2 = 2nah(2n + 1) 



(15) 



where the quantum number n is an odd integer (n — 
1, 3, 5, • • •) from the boundary condition: x(0) = 0. This 
result is identical to that obtained by Koma, Suganuma 
and Toki 0. 



B. 1^0 

Next we consider the case where the angular momen- 
tum does not vanish. Since the eigenvalue equation can- 
not be solved analytically, we make use of the WKB ap- 
proximation |13| . We put the wave function <p(r) as 



tp(r) = A(r)e lS(r ^ h 



(16) 



where an amplitude A(r) and the phase S(r) are un- 
known real functions. Substituting it into Eq.(13), we 
get 

H(r, e- iS / A pe iS / A , l)A{r) = EA{r) (17) 

Here we expand the left-hand side of this equation with 
respect to U. The lowest order equation leads to 



H(r,p,l,) = E, 



(18) 



where p stands for From this equation, we get the 
phase function S(r) and the energy eigenvalue E; This 
procedure will be done in Sec. 5. The next order of Ti gives 



2 ( ^ + ^ (r)= ° 



(19) 



where the operators on the left-hand side have been sym- 
metrized because it should be Hermitian. With the sim- 
ple calculation, this equation can be transformed into 

d ,dH 



dr^ dp ^ 



(20) 



Therefore the amplitude function A(r) is expressed by 

C 



A( r y 



dH r 2 
dp 



(21) 



where the C is a constant which should be determined 
by the normalization condition Eq.(10), that is, 



(22) 



Thus we have obtained the wave function of our system. 
The amplitude will be used to determine decay widths of 
glueballs, which is done in the next section. 



IV. DECAY WIDTHS OF GLUEBALLS 

As mentioned in the Introduction, the existence 
of glueballs are predicted theoretically and has been 
"found" on a lattice ^4|-[24|]. It is, however, not identi- 
fied experimentally. This may suggest that even if they 
are produced, they decay rapidly. Therefore it is very in- 
teresting to investigate the decay widths of our flux ring. 
To this end, we must introduce new interaction in our 
present model. 

We assume that this new mechanism is derived by the 
pair creation of qq: string — ► qq — > many hadrons. This 
mechanism that is called the Schwinger one 01 was used 
to calculate the decay widths of mesons and baryons in 
Rcf.(3) by us. The probability of the pair production in 
a unit space-time volume in the flux tube is given by 



2 00 1 



-), 



(23) 



q n—1 



where a is the string tention and indicates a sum- 
mation over all quark flavors with mass m q (q=u,d,s) 
19]- 2]J. The probability of the decay (pair production) 
in the time interval dt is given by 



dW = 



<p(r) \ 2 w£l(r)r 2 drdt, 



(24) 



where f2(r) denotes the volume of the ring when its radius 
is r. It is given by 



ft(r) = AS / y/l - v 2 rd(p 



(25) 



Substituting this equation into Eq.(24), we obtain the 
decay width, 



r = ^ = J™ Abr^l-{^fK 2 {k)A 2 dr (26) 
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where K 2 (k) is the complete elliptic integral of the second 
kind, which is defined by 

w / 2 . 

K 2 (k) = / J 1 - k 2 sin 2 <pd<p (27) 
Jo 



The new quantity k (0 < k < 1) is defined by 



k 2 = 



r 2 cv 2 



1 



(f) 2 



(28) 



The constant b — wASq defined in the right-hand side 
of Eq.(26) means probability of the string breaking per 
unit length per unit time. 

The decay width obtained above can be related to the 
excited energy in the case of zero angular momenta (I = 
0). First we notice K 2 (k = 0) = tt/2 due to u = 0. Then 
Eq.(26) is transformed into 



r = 



Air 2 ab 
E 



A 2 r 4 dr 



A-K 2 ab 
E 



(r 2 ) 



(29) 



The (r 2 ) denotes the expectation value for the excited 
state. Therefore the decay width is simply expressed by 



T = —E, 

2a ' 



(30) 



where E is the energy of the corresponding excited state. 
The decay width is proportional to the energy in the 
glueball spectroscopy. This fact is natural, because both 
the pair production and the energy are proportional to 
the volume of the flux tube. 



V. NUMERICAL RESULTS 

Since the mass spectra and the decay widths of glue- 
balls have been given in our model, we will carry out the 
numerical calculations, in particular in the case of the 
rotational motion. First we take up the energy and the 
angular momentum of our system, which are constants 
of motion, 



E = 



Aar r' 2 
1 ^-(W'o 



Aar 
Aar 3 u> 



\/l — k 2 sin 2 if 



dip 



it/2 



(31) 



/ 

J o 



sin 2 <p 



V~l - (W Jo v/l-^sinV 
Aar 3 uj K 1 (k)-K 2 (k) 



dtp 



k 2 



(32) 



where K\{k) is the complete elliptic integral of the first 
kind. If the first equation is squared, the energy is rep- 
resented as 



E z 



p 2 + \Qa 2 r 2 Ki{k) 2 = p +V 



(33) 




FIG. 2: The potential for the breathing mode. The L in the 
figure denotes for the angular momentum of the closed string. 



The V in the right-hand side is defined by 
V{r) =4arK 1 (k), 



(34) 



which means a scalar potential of the closed string. Here 
it should be noted that the potential is not at all linear 
function of r, because the k in this equation is a compli- 
cated function of r as shown in Eq.(32). 

It is instructive for us to figure out the potential V(r) 
as a function of r. From the other equation (32) on 
the angular momentum, the following relation can be de- 
rived, 



1 



Aak 
~T 

4a 



it/2 



sin 2 (p 



\/\ — k 2 sin 2 tp 



dtp 



= -( Kl (k) - K 2 (k)). 



(35) 
(36) 



This equation gives the function k(r) of r so that we are 
able to draw the graph of the potential V{r). The result 
is shown in Fig. 2. The repulsive core appears near the 
origin, which is due to the centrifugal potential. On the 
contrary, the potential becomes linear at larger distance. 
These behavior is natural and understood analytically as 
follows. Using Eqs.(34) and (36), the potential is rewrit- 
ten as 



hi 

V(r) = -+AaK 2 (k)r 
r 



(37) 

Then 

E -» 



When r — > 0, we note that k — > 1 from Eq.(35). 
the potential is reduced to the centrifugal one: 
yV + (l/r) 2 w p + l 2 /(2pr 2 ). On the other hand, it 
approsches the linear potential at larger distance, which 
is the confining potential: V(r) — * 2irar. 

The excited energy can be calculated using Eq.(33). 
From the boundary condition of the wave function, we 
obtain the Bohr-Sommerfeld formula: 

<j> p{r)dr = (2n+ 1)2^ (38) 
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FIG. 3: The mass spectrum of the vibrational mode as a 
function of the vibrational quantum number. 
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FIG. 4: The mass spectrum of the rotational mode as a func- 
tion of the angular momentum. 

If Eq.(33) is substituted into the left-hand side of this 
equation, the energy E is determined (n = 1, 3, 5, • ■ ■). 

There are two parameters in our model: the string 
tension a and the string breaking probability b. We take 
the same values as those used in the previous paper 
which reproduce the string tension a = 0.15 GeV 2 and 
the decay width T p = 151 MeV. 

First let us discuss the mass spectra of the excited 
states. The mass spectrum of the vibrational mode given 
from the Eq.(15) is shown in Fig. 3. The horizontal axis 
is the vibrational quantum number and the vertical one 
the squared mass. The quantum number must be odd 
integer (n = 1,3,5, •) as noted above. This straight line 
is an analog of the Regge trajectory of rotational mo- 
tion. The vibrational modes on the straight line cor- 
respond to the 1=0 states that lie on the (Regge) parent 
and (many) daughter trajectories with equal spacing. On 
the other hand, the spectrum of the rotational one cal- 
culated by Eq.(38) is shown in Fig. 4. The horizontal 
axis is the angular momentum (L — 0, 2, 4, • • •), In this 
case, the boundary condition at r = does not exist due 
to the existence of the centrifugal potential. The mass 
spectrum in Fig. 4 is drawn in the case of n — 0. We 
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FIG. 5: The decay widths of the vibrational mode as a func- 
tion of the vibrational quantum number. The triangles in the 
figure denote the ratio of the decay width to the mass. 
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FIG. 6: The decay widths of the rotational mode as a function 
of the angular momentum. The triangles in the figure denote 
the ratio of the decay width to the mass. 

have obtained an almost straight line: M 2 oc L. The 
Regge slope dL / dM 2 is smaller than that of the ordinary 
Regge trajectories of mesons and baryons. This discrep- 
ancy seems to originate from the geometrical structure 
of the string. If the total length of the string is fixed, 
the moment of inertia of the ring is smaller than that of 
the stick-like string (meson and baryons). On the other 
hand, the energy of the string would be proportional to 
the total length. Therefore the slope of the Regge trajec- 
tory of the glueball is smaller than that of the ordinary 
hadrons. This problem will be discussed later. 

The decay widths of the glueball are shown by the 
closed circles in Fig. 5 and Fig. 6. They are rather larger 
than those of the ordinary hadrons so that the glueballs 
generally decay easily. This fact explains that no glue- 
ball have been found experimentally. On the other hand, 
recent lattice QCD simulations suggest that the decay 
width of the scalar glueball is narrow [2^, I23L |24| , which 
is caused on the OZI rule. If the probability of the de- 
cay to one qq pair state is denoted by wi, we should use 
the w — w\ instead of w in our above calculations due 
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TABLE I: Mass spectrum of the blueballs. They are calcu- 
lated by the lattice QCD and our model. 





lattice QCD (MeV) 


Our model (MeV) 


Oi 


1.7 


1.68 


2 


2.5~ 2.8 


2.57 


2i 


2.4 


2.69 



state is the glueball state by the rotational motion and 
the energy is slightly smaller than that by our model. 
This is caused on the assumption that the shape of the 
ring is always circular. The centrifugal force would de- 
form the shape so that the ring becomes elliptic. Taking 
this effect into account, the moment of inertia of the ring 
becomes large and the excitation energy is lowered. 



to the OZI rule. Therefore the resulting decay rate is 
lowered. The triangles in the same figures represent the 
ratio of the decay width to the mass. In the case of the 
vibrational mode (Fig. 5), the decay widths are propor- 
tional to their mass. This is an exact result as shown in 
Eq.(30). The constancy of the ratio seems to hold also 
for the rotational mode as seen in the Fig. 6. 

The relation between the decay width and the mass is 
understood as follows. The mass of hadrons is propor- 
tional to the length of the string. On the other hand, the 
breaking probability is also proportional to the length of 
it because the cross section of the tube is the same for all 
hadrons (universality). Therefore we obtain the linearity: 
r cx E. 

Finally we add some comments on the experimental 
data [25|. There are some candidates: /o(1500) and 
/o(1710) for the scalar glueball; / 2 (2300) and / 2 (2340) 
for the tensor glueball. Unfortunately the abundance of 
q — q meson states in the 1-3 GeV region and the possi- 
bility of quarkonium-glueball mixing states still make it 
difficult to identify the glueball states. Instead we refer 
results by the lattice QCD simulation j(| , which is shown 
in Table I. The first and second J = states are expected 
to be the breathing modes of the glueball. These values 
are consistent with those by our model. The first J = 2 



VI. CONCLUDING REMARKS 

We have developped the flux tube model of glueballs in 
the context of that of mesons and baryons. The glueball 
is described by the closed string whose properties are the 
same as those of the stick-like string that constitutes a q— 
q meson or a q — (qq) baryon. Our model is characterized 
only by two parameters, the string tension a and the 
string breaking parameter b. If the two parameters are 
taken so as to reproduce the mass and the decay width of 
the rho meson, our model has no free parameter. Then we 
have investigated the mass spectra and the decay widths 
of glueballs with the use of the flux tube ring model. The 
numerical results are consistent with those of the lattice 
QCD. The calculated decay widths are sevral hundreds 
MeV. These values are rather large so that it may be 
difficult to identify the glueball experimentally. 
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